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Abstract Most statistical tools used to characterize the complex structures of the interstellar medium can be related to the 
power spectrum, and therefore to the Fourier amplitudes of the observed fields. To tap into the vast amount of information 
contained in the Fourier phases, one may consider the probability distribution function (PDF) of phase increments, and the 
related concepts of phase entropy and phase structure quantity. We use these ideas here with the purpose of assessing the ability 
of radio-interferometers to detect and recover this information. By comparing current arrays such as the VLA and Plateau de 
Bure to the future ALMA instrument, we show that the latter is definitely needed to achieve significant detection of phase 
structure, and that it will do so even in the presence of a fair amount of atmospheric phase fluctuations. We also show that 
ALMA will be able to recover the actual "amount" of phase structure in the noise-free case, if multiple configurations are used. 
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1. Introduction 

The physics of the interstellar medium (ISM) stands at the 
crossroads of many astrophysical problems, from stellar for- 
mation to galaxy evolution. Without a proper understanding of 
the processes taking place in the ISM, and of their interplay, 
complete and satisfactory solutions to these problems cannot 
hope to be met. 

Turbulence is one such process (se e 
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Consequently, a quantitative description of these struc- 
tures is a necessary first step towards understanding the 
physics of the ISM, and many statistical tools have 
been used to this e nd. Let u s me ntion th e powe r spec- 
trum (see e.g. iGautier et alJ [l992; Dickev et alJ 1200 it 
IStanimirovic & Lazman|| 200l|). the autocorr elation function 
/Kleiner & Dickman"l985';'Perault et al."l98^, the A-variance 
tStutzki et al. 1998; Bensch et al. 2Q01), the fractal dimen- 
sion Ipalgarone et alJ]l99ll) and the wavelet decomposition 
jCiill & Henriksenll99Ch . 

These various tools are not altogether independent from 
one another. By definition, the autocorrelation function is the 
Fourier transform of the power spectrum, to which the A- 
variance and fractal dime nsion can also be related, albeit 
less directly dStutzkiet alj ri998). Finally, the A-variance can 
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be written as the v ariance of wavelet transform coefficients 
dZielinskv & Stutz5itl999) . On the whole, it is then fair to say 
that all of these tools are connected, in one way or another, 
to the power spectrum, although some of them are of easier 
and more reliab le use depending on the type of observation 
(Stut zki etal J 19 98; Bensch et al. 2001). 

Since the power spectrum is given by the squared ampli- 
tudes of Fourier components, it basically ignores any struc- 
tural information that may be contained in the Fourier phases. 
Now, each Fourier component corresponds to a plane wave in 
direct space, with a given wave vector, amplitude and phase. 
The Fourier transform being linear, the observed structures 
are the result of the interaction between the various plane 
waves. Ignoring the phases when characterizing the structures 
is thus comparable to ignoring the interference phenomenon, 
and therefore marks a major loss in structural information. 
This has been confirmed by simple numerical experiments 
(Juvells et al. 1991; Coles 2005.) . 

In the experiment performed bv lCoIes| ( l2005l) . the Fourier 
phases of a numerical simulation of galaxy clustering, which 
is a highly-structured field, are randomly reshuffled in Fourier 
space. The resulting field has lost most of the filamentary struc- 
ture observed in the original image. This shows that it is in the 
Fourier-spatial distribution of the phases, and not in their values 
themselves, that most of the structural information must lie. 

The importance of this information may be best estimated 
in the context of interferometry. Indeed, interferometers es- 
sentially measure some Fourier components of the observed 
structures, and thus theoretically provide direct access to their 
phases. With the forthcoming ALMA instrument, the capacity 
of interferometers to detect structure in the Fourier phases in 
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real time may be assessed. This is the purpose of this paper, 
which is organized as follows: Section goffers a summary of 
the Fourier phase analysis technique, whose numerical imple- 
mentation is presented in section|31 The main part of the paper, 
dealing with the application of these techniques to interfero- 
metric observations, is the topic of section|3 Finally, section|5l 
gives a summary and conclusions. 

2. Fourier phase analysis 

The importance of Fourier phases in terms of structure 
has been exploited in various studies concerning vari- 
ations of the magnetic field in cometary plasmas and 
Solar wind (Polvaiannakis & Moussas 1995), the large-scale 
clustering properties of the Unive rse ( Sche rrer et 



large-scale 
t^ll99lt 



IChiang & CoiesI l200(i IWatts et alJ [2003^: £oles' '2005"), and 
Cosmic Microwave Background (CMB) maps (Coles etai] 
120041) . 

As we noted earlier, it is in the Fourier spatial distribu- 
tion of phases that information should be sought. To quantify 
this, IScherrer et alJ (Il99ll) suggested considering the statistics 
of phase increments As<p(k) = (p(k + 6) - <p{k) between points 
separated by a given lag vector S in Fourier space. 

In a field with no phase structure, all phases are uncor- 
rected, and phase increments are therefore uniformly dis- 
tributed random variables. Conversely, if the probability distri- 
bution function (PDF) of phase increments deviates from uni- 
formity, this deviation may be seen as a signature of phase 
structure. As an example, consider the column density of a 
5 12-' weakly compressibl e hydrodynamical turbulence simu- 
lation (> Porter etalJll994 . shown on Fig. ^ Since this field 
is periodic, its Fourier phases can be computed using a Fast 
Fourier Transform (FFT) algorithm. Given a lag vector 6, the 
PDF p {As<p) of phase increments for this lag is then approxi- 
mated numerically by computing the histogram of A^^ values. 
Two of these histograms are shown on Fig. |2l for lag vectors 
Cv and Cy, which are the unit vectors' of the Fourier space basis 
associated with a direct space basis (u^, m,.). We observe that 
these distributions are not uniform, with a single wavelike os- 
cillation around the value 1 /(2n), and that the amplitude of the 
oscillation is more important for 6 = e ^ than for 6 = e, , which 
may be interpreted as evidence for anisotropy. 

For other lag vectors or diff'erent images, the shape of 
the histogram remains, although the amplitude may vary. As 
pointed out bv lWatts et al.l ll2003( ). the underlying distributions 
are very likely to be von Mises distributions, given by 



1 



exp [-KS cos (Ag^ - fig)]. 



where /o is the zeroth order modified Bessel function of the 
first kind, and the parameters ns and ks > control respectively 
the position of the distribution's minimum and the oscillation's 
amplitude. Thus, Kg can be viewed as a measure of the amount 
of phase structure in the image. For the histograms of Fig.|2l 



' Their lengths are actually l/A'^-A, and l/A'^A,,, where A,, x A,, is 
the actual size of a pixel in direct space and A',- x A',, is the image size 
in pixels. 




X (pixels) 

Figure 1. Column density of a 512^ weakly compressible hy- 
drodynamical turbulence simulation obtained by Porter et alJ 
(Il994l) . used here as a model brightness distribution for phase 
structure analysis. 



on which fits by von Mises distributions are shown, the values 
of Ks found are /fg, = 0.197 and Ke^^ - 0.0656. 

Phase entropy, introduced by 'Polvgia nnakis & M oussa^ 
( 1 995 ) . is another measure of the distribution's departure from 
uniformity, and is defined by the integral 



S{6) 



-f 

J -n 



p(A50)ln[p(A50)]dA5 



It reaches its maximum value Sq - In (2:7r) for uniform PDFs, 
and tends to -oo for 5-function PDFs. These extreme cases 
correspond respectively to fields with purely random phases 
such as fractional Brownian motions^ (see e.g. Istutzki et all 
Il998l) . and to fields containing a single point-source. Given 
these limits, it is convenient to consider the positive quantity 
Q{6) - Sq - S(6), which we dub "phase structure quantity" 
in the rest of this paper There is a monotonous relationship 
between the two measures a-^ and Q(6), given by^ 

but we shall deal only with the phase structure quantity Q from 
now on. The reason for this is computational: phase struc- 
ture quantities can be computed on histogram-like functions, 
while the parameter a:^ requires fitting them by von Mises 
distributions first, a procedure which may not converge prop- 
erly when the underlying distribution is close to uniform. For 
the von Mises distributions fitted on the histograms of Fig.|2j 
<3(c,) = 9.6 10"^ and Q{ey) =1.1 10 ^ respectively. Using col- 
umn density fields of the same compressible turbulence sim- 
ulation at different times, we find that the typical values for 



- These are random fields characterized by a power-law power spec- 
trum and random phases. 
^ /i is the first order modified Bessel function of the first kind. 
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Figure 2. Histograms of phase increments for the turbulent 
brightness distribution (Fig.[3. The lag vectors used are 6 - e ^ 
(top panel) and 6 - e,. (bottom panel). The dotted lines repre- 
sent the uniform distribution and the dashed lines are fits by von 
Mises distributions (see text). The number of bins is n = 50. 



phase structure quanti ties lie typically below 10 ^, whereas 
IChiang&Cole3ll200(t found values as high as ~ 0.4 for grav- 
itational evolution of density perturbations. Although density 
contrasts are quite different, making direct comparison some- 
what hazardous, this suggests that Fourier phase analysis may 
be useful in determining the physical processes governing the 
formation of structures in the interstellar medium. 

3. Phase structure quantity in practice 

Due to the limited number of phase increments that can be 
computed from a finite-sized image, the histograms do not per- 
fectly sample the underlying PDFs, but rather include a cer- 
tain amount of statistical noise. Consequently, the histograms 
of phase increments for images such as fractional Brownian 
motions are not exactly uniform, and the corresponding phase 
structure quantities are not zero. More generally, the phase 
structure quantities Q associated with the underlying distribu- 
tions are to be distinguished from those found by numerical in- 
tegration performed on the histograms, since these depend on 
the number p of increments available and the number n of bins 



used. We write Q for these estimates of phase structure quanti- 
ties. For the histograms shown on Fig.|2] built with « = 50 and 
/7 = 512 X 511, we find(3(e,) = 9.8 10"^ and(3(ev) = 1-3 10-^ 
which shows that the difference with Q (about 2.10"'*) can be- 
come significant for low phase structure quantities. 

Assessing the detectability of phase structure makes it nec- 
essary to determine the threshold of Q above which it can be 
said that an image deviates significantly from a "structureless" 
field, given values for n and p. Mathematically speaking, this 
amounts to determining an upper limit to the probability that 
the estimate Q of the phase structure quantity be greater than a 
given positive real number x, assuming a uniform parent distri- 
bution, as a function of n and p. 

We may obtain such an upper limit by an analytic approach, 
as described in detail in t he appe ndix. The demonstration is 
based on results obtained bv lCastellan C2000) in her PhD thesis, 
which is avail able onUne"*. These resu lts are themselves derived 
from those of iBarron & Sheul ( Il99lh and show that an upper 
limit to the probability P > is given by the quantity 



1 - Erf e 



2(n - 1) 



2 ^ 2px{\ - ef 



l+e 



where Erf is the error function and;^'^ is the chi-square statistics 
of degree n - 1. The first term Pi corresponds to the probability 
of having an histogram with a large fluctuation, the meaning 
of "large" being defined by means of the arbitrary positive real 
number e, as explained in the appendix. The second term P2 
corresponds to the probability of having Q > x with a more 
regular histogram. Since the threshold value x appears in P2 
only, we should look for a value of e ensuring that Pi <s; P2, so 
that we may decide on the phase structure quantity threshold 
using well-known chi-square statistics. 

This suggests an "e-adaptive" procedure, which is the fol- 
lowing: For each {n,p) pair, the value of e is chosen so that 
Pi is small. The quantiles of the chi-square statistics are then 
used to extract a value of x so that Pi <«c P2 <sc 1. In practice, 
we took Pi = 10 and P2 = 10"^. In the end, there is a 0.99 
probability that the phase structure quantity be less than x with 
these values of n and p, assuming that the underlying distri- 
bution is uniform. Conversely, if the measured phase structure 
quantity Q is greater than x with these values of n and p, then 
the underlying distribution is most likely non-uniform. 

As will be clear later on, this procedure may lead to very 
conservative upper limits. That is why we may also follow a 
"fixed e" procedure, in which we set e to a fixed value, say 0.1, 
and simply ignore Pi. As in the e-adaptive procedure, we then 
compute X so that P2 — 10"^. 

The influence of the number of phase increments p and the 
number of bins n on the reliability of the phase structure quan- 
tity may also be studied numerically. To this end, we have com- 
puted a series of two-dimensional fractional Brownian motions 
of various sizes. The lag vector being fixed, namely S - e^, 
the numerical calculation of the phase structure quantity for 
these simulations provides us with an estimate of the limit 
above which one should conclude that phase structure is in- 
deed present in the image. More precisely, we have computed 

^ _http: //www .math. u-psud.fr/theses-orsay/2SSS/6Q39 .html I 
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Figures. Evolution of phase structure quantities Qie^) for 
fractional Brownian motion fields of size pi x pi (implying 
p = Piipi - !)■ The symbols represent the mean values of the 
phase structure quantities for ten realizations of the fractional 
Brownian fields, while the vertical lines represent the standard 
deviations. The different symbols correspond to n - 50 (dia- 
monds) and n - 500 (squares). The dotted and dashed lines 
represent the corresponding theoretical upper limits, according 
to the e-adaptive procedure (grey lines) and the fixed e proce- 
dure (black lines). The cross represents Q(e^) in the case of the 
turbulent column density of Fig.[n for n - 50. 



fractional Brownian motion fields of size varying from 32 x 32 
to 800 X 800 pixels^. For each size, ten fields were built and 
two histograms drawn from the maps of their phase increments, 
with respectively n = 50 and n = 500. For each pair (n, p), 
we then computed the mean and standard deviation of the ten 
phase structure quantities associated with the fields. The results 
are shown on Fig.|3 

It appears that the computed phase structure quantity in- 
creases as the size of the image decreases, and as the number 
of bins increases. This is interpreted by the fact that the his- 
tograms are then less accurate samples of the underlying dis- 
tributions. The figure also shows theoretical upper limits com- 
puted using both procedures described earlier. Unsurprisingly, 
the e-adaptive upper limits fall above the values found in the 
numerical simulations, even well above them, which demon- 
strates the conservativeness of this approach. On the contrary, 
fixed e upper limits for n = 50 match the numerical esti- 
mates better, but actually fail when n - 500. This is due to 
the fact that, in this latter case. Pi > 1, which makes the upper 
limit thus computed useless. The position of the phase structure 
quantity Q{ex) for the turbulent column density (Fig.^ in this 
plot is quantitative evidence that this field does harbour phase 
structure. 



4. Application to interferometric observations 

4.1. Introduction 

In the ideal case, interferometers sample the Fourier trans- 
form of observed brightness distributions, and therefore al- 
low direct measurement of the phases of Fourier components. 
Consequently, it is theoretically possible to have access to 
phase increments, and to phase structure quantities for the 
observed fields. Since this can be done in real time, as the 
Earth's rotation allows for a better sampling of Fourier space 
( [Thompson et al.ll99ll) . we may look for the minimum observ- 
ing time required to detect a significant phase structure quantity 
in the data. This is the topic of this section^] 

Before we carry on, however, we should stress that, in prac- 
tice, phases measured by the instrument do not directly yield 
the actual phases of the model brightness distribution. First 
of all, the antennae are not pointlike receptors, so that bright- 
ness distributions are multiplied by a primary beam. In Fourier 
space, this corresponds to a convolution of the Fourier com- 
ponents by a finite size kernel. Consequently, the phase mea- 
sured at a given point in the (m, v) plane does not yield the ac- 
tual phase of the Fourier component of the brightness distribu- 
tion at that point, but involves all Fourier components within 
a small neighbourhood. One may circumvent this difficulty by 
considering mosaicing observation techniques. In short, these 
amount to imaging large fields by pointing the array towards 
differe nt directions succ essively, and "gluing" the subfields to- 
gether dBhatnagar et al.ll?005ii) . In image space, this last step is 
done in such a way that the fall-off due to the primary beam 
pattern in a subfield is compensated by the rise of the primary 
beam pattern in the adjacent subfield, so that the effective pri- 
mary beam is more or less uniform over the large composite 
field. In Fourier space, these techniques correspond to a finer 
sampling of the (m, v) plane, effectively reducing the size of the 
convolution kernel, so that the measured phases are better es- 
timates of the actual phases of Fourier components. We will 
not discuss this problem any further here, as it should require 
an extensive study that is not within the scope of this paper, 
and we assume from now on that antennae can be modelled as 
pointlike receptors. 

Another point to consider is the fact that measurements in 
the (m, v) plane are not sampled on a regular grid, and regrid- 
ding is widel y used to allow for Fast Fourier Transforms to 
be performed dThompson & BracewelL1974.) . This means that 
the measured phases end up being associated with a different 
wavenumber than the one they actually correspond to. In our 
case, as we use model brightness distributions that are already 
sampled on a regular grid, the gridding problem can be by- 
passed^. 

Lastly, noise contributions, especially those due to the tur- 
bulent fluctuations of the atmosphere, blur the true phase val- 
ues. This effect will be discussed in detail, in section |431 



^ So that p varies from 32 x 31 to 800 x 799. 



^ This is actually only valid in the noise-free case (see l4.5l and l4.6> . 
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Table 1. Instrumental parameters 



Wavelength 


Longitude 


Latitude 


Dump time 


1.3 mm 


-67.75 - 


-23.02 ° 


10 s 



Table 2. Array characteristics. For each configuration, the min- 
imum and maximum separation between any two antennae are 
given. 



Instrument 


ALMA 


PdB 


VLA 


Antenna diameter (m) 


12 


15 


25 


Number of antennae 


60 


6 


27 


A configuration (m) 


19- 11527 


32 - 400 


807 - 37235 


B configuration (m) 


76 - 3005 


71 - 331 


247 - 11314 


C configuration (m) 


83 - 2303 


48 - 229 


79 - 3444 


D configuration (m) 


43-1618 


24-113 


41 - 1048 


E configuration (m) 


34 - 909 


n.a. 


n.a. 


F configuration (m) 


15 - 229 


n.a. 


n.a. 




200 300 400 

X (pixels) 



500 



4.2. Simulations of observations 

The instrument simulator used is of the simplest kind, and its 
parameters, taken to match those of ALMA, are summarized 
in Tabled The instrument tracks the source as long as it re- 
mains above a minimum elevation^ of 10°, which, given the 
array's latitude and the source's chosen declination of -20°, 
represents a maximum integration time of 11 hours and 38 
minutes. Regarding the number and positions of the antennae 
on the groun d, we have chosen configurations optimized by 
iBoondOOOll) based on ALMA specifications. For comparison, 
we have also considered configurations taken from current ar- 
rays, such as the Plateau de Bure (PdB) radio-interferometer 
and the VLA. To make meaningful comparisons, we have used 
fictitious arrays located at the same geographical coordinates 
as ALMA, observing the same source. Only the number and 
positions of the antennae are changed to match the configura- 
tions of the PdBI and the VLA. The characteristics of all the 
arrays used are summarized in Table|2| Using the source's ap- 
parent movement in the sky and the locations of the antennae 
on the ground, we obtain ungridded (m, v) covers as functions 
of integration time. 

As model brightness distribution, we use the turbulent col- 
umn density shown on Fig.^ which we know harbours phase 
structure (Fig. O. For comparison purposes, a field with the 
same power spectrum and random phases is also considered 
(Fig.l^. Both model fields are 512x512 images, so (u, v) cov- 
ers are regridded on a grid of that size, using nearest-neighbour 
interpolation. The size of the (m, v) cells is chosen to be half the 
antenna diameter, to satisfy the Nyquist criterion^. This limits 
the size of the maximum baseline that can be considered. Only 
the F and E configurations of the ALMA instrument, the D con- 
figuration of the VLA instrument and all configurations of the 

' No shadowing of the antennae is taken into account. 
This means that the actual pixel size is dilferent for the three arrays 
considered. We shall not be troubled by this, given that the brightness 
distributions used are scale-invariant and are not subject to boundary 
conditions, so that their actual physical extents need not be specified 
and can therefore be scaled accordingly. 



Figure 4. Synthetic field with power spectrum identical to that 
of the turbulent column density field (Fig.QJ, and fully random 
phases. 

PdB instrument fiton512x512 grids with the corresponding 
pixel sizes. These are used in the single-configuration simu- 
lations. However, it is possible to consider the more extended 
configurations, provided the longest baselines are ignored. We 
did so in the case of multi-configuration observations, as ex- 
plained in section l4!4l As an illustration, Fig.|5lshows the grid- 
ded covers for the F and E configurations of the ALMA instru- 
ment, considering the time period of one hour centered on the 
source's transit at the meridian. 

4.3. Evolution of the measured phase structure with 
integration time 

As the observation is carried out, more and more Fourier phases 
are measured and the number p of phase increments increases, 
for any given lag vector 6. The question is whether this allows 
to bring the theoretical and numerical upper limits of Fig. |3] 
down sufficiently, below the measured phase structure quan- 
tities, to ensure positive detection. To answer this question, 
Figures|6lto|9lshow the evolution of the measured phase struc- 
ture quantities as a function of integration time. 

Concerning the E configuration of the ALMA instrument 
(Fig. |6j, the conclusions that can be drawn are the following: 
For 6 - e I, a short integration time of approximately twenty 
minutes is enough to conclude that phase structure is present 
in the image, since the measured value then becomes larger 
than the adaptive procedure's upper limit. On the other hand, 
the phase structure for 6 - Cy is harder to extract, due to the 
lower value of Q(ey) for the complete field. Approximately 7.5 
hours of integration are required to see the measured phase 
structure quantity rise above the more conservative theoretical 
upper limit, although the curves for the turbulent field and its 
random phase version are clearly distinguished on the whole 
range plotted. When comparing the two panels of Fig.|6l it ap- 
pears that the measured phase structure quantity for the max- 
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Figure 5. Gridded (u, v) covers for the F (upper panel) and E 
(lower panel) configurations of the ALMA instrument, when 
observing the source during one hour centered on the meridian 
transit. Pixel sizes are the same in both cases, namely 6x6 
meters, and the F configuration cover has been zoomed. The 
maximum pixel values are 1726 samples per cell in the F con- 
figuration and 470 samples per cell in the E configuration. 



imum integration time is larger than the value for the whole 
field in the case 6 = c v, while it is the opposite for 6 = e, - 
These discrepancies are due to the limited range of measured 
spatial frequencies, and are addressed in section We point 
out, however, that the approximate constancy of the measured 
phase structure quantities for T > 3.10"* s shows that the (m, v) 
plane is then close to being fully sampled within this range. 

Regarding the F configuration of the ALMA instrument 
(Fig.0, its compacity leads to a smaller number of phase in- 
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Figure 6. Evolution of measured phase structure quantities 
with integration time, for the E configuration of the ALMA 
instrument. The upper panel corresponds to 6 - Cx and the 
lower panel to 6 = Cy. The black solid lines correspond to the 
turbulent model brightness distribution (Fig.Q, while the grey 
solid lines correspond to the random-phase brightness distri- 
bution (Fig.E). The dotted lines represent the phase structure 
quantities of the complete turbulent brightness distribution, for 
each lag vector, and the dashed lines represent the evolution of 
theoretical upper limits with integration time, using the fixed 
procedure (black lines) and the adaptive procedure (grey lines). 



crements, which makes phase structure detection all the more 
difficult. Our best chances lie with the S - lag vector (up- 
per panel). In this case, the measured Q remains constantly be- 
low the upper limit set by the adaptive procedure, but becomes 
larger than that set by the fixed e procedure, for an integra- 
tion time of just above one hour. This may be seen as evidence 
that the phase structure can also be detected with this compact 
configuration. Whether it should be possible to determine the 
phase structure quantity for the whole field is uncertain,since 
it falls below the less conservative upper limit, but above the 
curve corresponding to a random-phase field. 

Figures |8] and |9] show the evolution of the measured phase 
structure quantity (3(e i) as a function of integration time, us- 
ing the B configuration of the Plateau de Bure interferometer 
and the D configuration of the VLA, respectively. What ap- 
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Figure?. Same as Fig. |6lbut for the F configuration of the 
ALMA instrument. The adaptive procedure's upper Umits lie 
outside of the range of Q values plotted. 

pears clearly on Fig.|8lis that the number of phase increments 
measured by the Plateau de Bure is insufficient to detect phase 
structure, as the curves for turbulent and random-phase bright- 
ness distributions are indistinguishable from one another. The 
same conclusion prevails for other configurations of this in- 
strument. On the contrary, the VLA allows such a detection, 
although it takes a long integration (about 6 hours) to see the 
phase structure quantity measured emerge from the adaptive 
procedure's upper limit. Let us note however that the curves 
for both model brightness distributions start going apart after 
less than twenty minutes of integration, which should give ob- 
servers a first hint that phase structure is present in the field. 
This diagnosis can be performed in real time by drawing ran- 
dom phases for the visibilities as they are measured. 

4.4. Multi-configuration observations 

Comparing Figs.|6lto|9l the "best" situation appears to be the 
E configuration of the ALMA interferometer, for 6 - e^. Yet, 
the phase structure quantity obtained after one transit of the 
source is not the one computed on the whole field. This is due 
to the fact that only 24% of the 512 x 512 Fourier phases are 
measured by this configuration. Using more extended configu- 
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Figure 8. Same as the top panel of Fig. |6j but for the B config- 
uration of the Plateau de Bure instrument. As with Fig.0 the 
adaptive procedure's upper limit lies outside of the range of Q 
values plotted. 
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Figure 9. Same as the top panel of Fig.|6j but for the D config- 
uration of the VLA instrument. 

rations, one should be able to recover the Fourier components 
lying outside the radius covered by the E configuration, and 
therefore hope to recover the correct value of the phase struc- 
ture quantity by combining visibilities from multiple configu- 
rations. Obviously, this can only be done after the observations 
have been carried out, unlike what has been considered until 
now, so that we should rephrase our initial question: what is 
the minimum integration time necessary in each configuration 
to ensure that combining the observed visibilities will lead to a 
"full" coverage of the (m, v) plane?' In this multi-configuration 
approach, we have to consider the extended configurations that 
do not fit on a 512 x 512 grid with the given pixel size, and 
simply ignore the visibilities falling outside the grid. 

Figure[roishows the evolution of the measured phase struc- 
ture quantity Q{ex) with integration time, using this approach 

' We should stress that we do not consider a set of integration times 
(ti , . . . , Tjv), one for each of the N configurations, but a single value t 
that applies to each configuration separately. 
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Figure 10. Evolution of measured phase structure quantity 
Q{ex) with integration time for an observation using all con- 
figurations of the instrument in turn. The black solid Une corre- 
sponds to the six configurations of the ALMA instrument, and 
the grey solid line to the four configurations of the VLA in- 
strument. The dotted line represents the value of v) for the 
whole field. 

with the ALMA and VLA instruments. In the former case, it 
appears that the phase structure quantity for the whole field is 
recovered with an integration time of 9 hours in each configu- 
ration, totalling 54 hours of observing time. On the other hand, 
a full-day integration using the four configurations of the VLA 
instrument is not sufficient to reach the phase structure quantity 
of the whole field. This emphasizes the better coverage of the 
Fourier plane that will be achieved by ALMA. 

4.5. Atmospheric phase noise 

In the previous subsections, it was assumed that interferom- 
eters sample the true Fourier phases of the observed fields. 
However, it is well known that turbulent motions within the 
atmosphere above the instrument alter the measured phases. 
Indeed, these motions cause the amount of water vapor along 
the line of sight to vary both in time and from one antenna to 
the other. This results in a delay error and thus a phase error 
for each baseline. This problem has been addressed thoroughly 
bv lLavl lll997albl) . and simulations of the atmospheric fluctua- 
tions at th e ALMA site of Chajnantor have been performed by 
LStirlin^ et al. ( 2005i) . Here, we have chosen to perform a sim- 
ulation of the effects of atmospheric phase noise by introduc- 
ing an atmospheric mask giving the refractivity field y, z) 
above the instrument. In practice, we assumed that this field 
can be regarded as a 200-m thick layer of frozen Kolmogorov 
turbulence, that is being transported along the east-west direc- 
tion at a wind speed of 2 m.s"' . Using a spatial resolution of 10 
meters, we computed as a fractional Brownian motion of size 
8500 X 150 X 20 pixels, with spectral index -1 1 /3. Integration 
of § along the different lines of sight for each antenna as the 
observation is performed yields phase delays, which are subse- 
quently correlated in order to give the atmospheric phase noise 
<pa{a,l3,t) for each pair of antennae (a,/!), at all times t. The 



field 1? is normalized so that the rms phase noise ctq for a pair 
of antennae observing the zenith and separated by a baseline 
d - 100 m should be one of a few specific values, namely 15°, 
45° and 90°. 

At each time step, a number of visibilities fall in each cell 
of the gridded (m, v) plane, and the observed "mean visibility" 
in the cell Ck, centered on wavevector k, is given by 



V'{k,t) 



1 



Nik, t) 



2 ^ V\k)cxp[i<f>,(a,pj')], (1) 



r'<t (afi/)fiCk 



where N{k, t) is the cumulative number of visibilities within 
the cell Ck at time f, and is the true visibility in this cell. 
This relation allows for the simulation of the measurement of 
Fourier phases in the presence of a given amount of atmo- 
spheric phase noise. Fig. II II shows the evolution of the mea- 
sured phase structure quantities Q{e^ and QiCy) in the case of 
the E configuration of the ALMA instrument. 
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Figure 11. Evolution of the measured phase structure quantity 
with integration time, in the presence of atmospheric phase 
noise (solid lines, with ctq specified next to each curve). The 
array used is the E configuration of the ALMA instrument, and 
the top and bottom panels correspond respectively Xo S - 
and 6 = c, . The dotted lines represent the phase structure quan- 
tities for the whole field, and the dashed lines correspond to the 
e-adaptive upper limits. 
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Provided that the value of the phase structure quantity is 
large enough for the whole field, as is the case for 6 = e^, 
the presence of phase structure can be easily detected in the 
presence of a fair amount of atmospheric phase noise. Indeed, 
even a rms phase fluctuation of 90° at 100 m is insufficient to 
bring the measured phase structure quantity down below the e- 
adaptive upper limit. In the case of 6 - e,, the conclusion is 
not so clear-cut: although the ctq - 15° curve exhibits values 
larger than the upper limit, it does not remain above this limit 
after a certain integration time. Faced with such a situation, one 
should be suspicious of the presence of phase structure in the 
ob served field'" . 

iButler et'aP (1200 ih performed measurements of the atmo- 
spheric phase fluctuations above the Chajnantor site, using a 
two-element interferometer observing a 11.198 GHz beacon 
broadcast by a geostationary satellit e positioned 35° above the 
horizo n. Using the scaling relation (lLavlll997al [Stirling et alJ 
l2005l) 

giving the rms phase delay as a function of the baseline length 
d, the wavelength A and the elevation angle their results 
translate to a noise level showing diurnal as well as seasonal 
variations going from cry ~ 14° to ctq ~ 57°. Consequently, 
the phase structure quantity for 6 = e ^ will undoubtedly be de- 
tected without any phase correction, although the use of dedi- 
cated water vapor radiometers, as is planned for ALMA, should 
allow for an effective decrease of the atmospheric phase noise 
by a substantial factor (Lav 1997b), making possible the mea- 
surement of the actual phase structure quantity for the whole 
field, using multiple configurations. 

4.6. Extensions 

It should be stressed that the approach used here is not to be 
considered optimal, but only as evidence that ALMA will be 
able to detect phase structure in the presence of atmospheric 
phase noise. Indeed, consider two tracks in the (m, v) plane, cor- 
responding to two different antenna pairs, going through the 
same grid cell. In the above approach, the measured phases 
in this cell may be very different for one antenna pair and the 
other, since these sample different lines of sight through the at- 
mospheric mask. As a result, the estimation of phase through 
the "averaged visibility" of Eq. (^ is badly contaminated. A 
more elaborate use of phase information would therefore have 
to be baseline-based. Keeping track of the phase measured by 
each baseline as a function of time, and computing phase incre- 
ments along the baseline's track should markedly reduce con- 
tamination by atmospheric phase noise. An important point in 
this approach, which is currently under study, is that the lag 
vector 6 is no longer a control parameter, but a function of time 
and of the baseline. 

In any case, phase increments corresponding to various lag vectors 
6 can be simultaneously computed, in real time, so that for a given 
observed field, it is enough that there exist one such lag vector for 
which phase structure detection is feasible. 



Conversely, the consideration of the radial evolution of 
phase structure quantities leads to another possible extension 
of this work, which is the inclusion of the kinematic dimen- 
sion. Velocity information is indeed accessible with high spec- 
tral resolution receivers such as those that will be used for 
ALMA (4096 spectral channels over 16 GHz bandwidths). 
Consequently, Fourier phase analysis applied to individual 
channel maps may prove a valuable tool for assessing the three- 
dimensional structure of velocity fields. We may for instance 
wish to compare phase structure quantities found across line 
profiles, and see if values found for individual channels are 
greater than those found for integrated emission maps, which 
is likely to be the case, as they present higher contrasts. 

5. Summary and conclusions 

In this paper, we have addressed the ability of radio- 
interferometers to detect and recover the information contained 
in the Fourier-spatial distribution of phases, which was previ- 
ously shown to store a vast amount of information about the 
structure of images. The PDF of phase increments and the re- 
lated concept of phase entropy were introduced in this per- 
spective. We ourselves have used the phase structure quantity 
Q, which is a minor modification of phase entropy leading to 
Q - for fields with purely random phases. 

Our main conclusion is that the dynamical range of spatial 
frequencies observed by the instrument is the key parameter 
allowing detection and measurement of phase structure. 

Using a turbulent model brightness distribution and instru- 
mental configurations based on the characteristics of the future 
ALMA interferometer and of two existing arrays (VLA and 
Plateau de Bure), we have assessed the minimum integration 
time required by each configuration to have a significant detec- 
tion of phase structure in the observed field. In the most conser- 
vative assessment, it appears that for a whole-field phase struc- 
ture quantity Q ^ 10"^, detection is achieved with a twenty 
minute integration in the ALMA E configuration (baselines go- 
ing from 34 m to 909 m), or with a six hour integration in the 
VLA D configuration, but is not achieved by any other instru- 
mental configuration tested". With a whole-field phase struc- 
ture quantity Q ^ 10"^, certain detection can only be achieved 
using the ALMA E configuration, in which case it takes about 
7.5 hours of integration. 

However, less conservative criteria allow for early hints at 
the presence of phase structure in the observed field. Indeed, 
by drawing random phases for the visibilities in real time, it is 
possible to compare the evolution of the phase structure quan- 
tity for the observed field to that for a random-phase field, and 
check if they start going apart at some point. This is the case 
for all ALMA and VLA configurations, with whole-field phase 
structure quantities Q ^ 10"^ and Q ^ 10"^, but not for any of 
the Plateau de Bure configurations. 

Regarding the possibility to recover the actual values of 
the phase structure quantity for the complete field, only multi- 
configuration observations with the ALMA instrument seem to 

' ' It should be reminded that the more extended configurations have 
not been used in this single-configuration approach. 
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allow for it, and it takes 9 hours in each of the 6 configurations 
to achieve this. 

Finally, we have studied the influence of atmospheric phase 
noise on the single-configuration observations, using the E con- 
figuration of ALMA and whole-field Q ^ 10"^. The maximum 
rms phase fluctuations that can be allowed without completely 
washing out the actual phase structure lie well above the typical 
range of variations for the Chajnantor site. Consequently, the 
use of water vapor radiometers to correct for the atmospheric 
phase fluctuations does not appear as a necessary feature of the 
ALMA array in this respect, although it should allow for a more 
accurate determination of the actual phase structure quantity in 
the multiconfiguration scheme. 

Possible extensions to this work include the study of phase 
increments along the baseline tracks, which should consider- 
ably reduce the effects of atmospheric phase noise, and the 
evolution of phase structure with frequency in high spectral- 
resolution observations of line sources. 
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Appendix A: Theoretical upper limits 

Independently of the underlying distribution, the parameters 
which have an influence on the histograms of phase increments 
are the number of samples p and the number of bins n. Assume 
then that p phase increments are drawn from a uniform distri- 
bution on ^ = L-7r,7r] and distributed over n intervals of 
equal length, thus yielding a n-binned histogram h. The phase 
structure quantity Q associated with h is nonzero due to the 
samphng noise. We therefore wish to obtain an upper Umit to 
the probabiUty P ({Q > x]j, as a function of x. The histogram 
values [hi] are normalized according to 



1=1 



Let us then define the functions 

1 " 
3(x) = — lyi and six) = Y^hilji„ 



i=l 



where the symbol l[a,h] stands for the function which is equal to 
one on the interval [a, b] and zero outside. The phase structure 
quantity associated with h may then be written as 

soix)ln[soix)]dx + I i(x) In [i(x)]dx. 

Now, since so is actually a constant on \-n,7t\, and since both 
s and sq are normalized to unity, it is straightforward to obtain 



Q 



£ 



s{x) In 



six) 



so(x) 



dx. 



which shows that the phase structure quantity may be inter- 
preted as a KiiUback pseudo-distance of * to io- The method of 



F. Levrier et al.: Fourier phase analysis in radio-interferometry 



11 



ICastellaiil ( I2OOOI) suggests to find an upper limit to P ({(3 > xfj 
by treating separately the cases of regular and extraordinary 
histograms, the latter being when the histogram presents an un- 
usually large or unusually low value. Let us then define, for any 
e > 0, the event 

In 1 
= {3/; - r| > er) with x,- = — /i; and r--. 

n n 

The event Q.^ is precisely the occurrence of an extraordinary 
histogram, and it depends on the real number e. For instance, 
if e = 0.1, this event occurs if one of the histogram values de- 
viates from the uniform value 1 l{2n) by more than ten percent. 
The value of x, may here be interpreted as the mean number of 
successful events in a series of p Bernoulli trials, the event in 
question, whose probability is r, being that a phase increment 
belongs to interval yi,. The usual values of p and n (in our case 
p - 512^ and n - 50) make it reasonable to assume that the 
central limit theorem applies. We may then write 

P(a) < nP ({\x, -r\> er}) ^ « [1 - Erf(x,)] , 

introducing the value of the error function 



Erf(;c) = ^ I e-''dt at 



Jo V 2(n - 1) 

Numerically, for p - 512^ and n - 50, the probability that 
there exist an interval yi, containing a number of increments 
different from the theoretical value by over five percent (e = 
5 10"^) is less than about 0.012, and it falls below 10"" for 
a ten percent discrepancy. In this case, it will be possible to 
neglect the contribution of extraordinary histograms. 

Regarding the regular histograms, for which does not 
occur, we have the following result, due to Castellan (2000) . 



mf(s. So) In I — j d.x < 2(3 < J sup(i, sq) In | 



2 

dx. 



For the regular histograms considered here, we obviously have 
(1 - e)so < inf(5, so) < sup(s, sq) < (1 + £)^^o, so that, using 



(1 + e)2 p 



I So In — 



1 

(1 - e)^ p 



which were obtained bv lCastellan ( 2000), we conclude that the 
phase structure quantity is equivalent to the statistics of de- 
gree'^ « - 1, for regular histograms, 

y^-y^ 1-e 1-i-e 
flf — < (3 < bf^ — with flf T- and b^ = 



2(1-1-6)2 ' 2(1 -e)2' 

For regular histograms, the probability we are concerned with 
is therefore subject to the inequality 



p({Q>4)<p({^.y >x|j = p(|;r 



2 2(1 - 6)^px ] 
1+6 I 



and taking into account both regular and extraordinary his- 
tograms, we come up with the upper limit given in the main 
body of the paper. 



The degree is n- 1 and not n because of the constraint that s should 
be normalized to unity. 



